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(54) METHOD AND APPARATUS FOR ELECTROCHEMICAL MEASUREMENT USING 
STATISTICAL TECHNIQUE 



(57) A predetermined voltage is applied to a biosen- 
sor twice to promote an electrochemical reaction, and 
the following parameters (P1 and P2) are calculated 
from the values of detected current. A statistical tech- 
nique is used with these parameters to compensate for 
errors so that the concentration of an object can be 
determined. P1: the ratio (If/lb) of the maximum current 
(if) or a current occurring the maximum in the first exci- 
tation to a current (lb) read at any point in the second 
excitation. P2: a current (lb) read at any point in the sec- 
ond excitation. 
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Description 
TECHNICAL FIELD 

5 [0001] This invention relates to a method for electrochemically measuring a concentration of an analyte in a sam- 
ple. More specifically, this invention relates to a method for electrochemical measurement to obtain a concentration of 
analyte in a sample by correcting errors caused by sample physical properties and sensor sensitivity. 

BACKGROUND ART 

10 

[0002] Recently in the fields of biochemistry, clinical medicine or the like, electrochemical measurement has been 
used to measure a concentration of an analyte in a sample in a rapid and simple manner. An electrochemical measure- 
ment is a method to measure an analyte in a sample by using an electric signal such as a current that is obtained from 
a chemical reaction or an enzymatic reaction. 

15 [0003] For example, the concentration of glucose in a blood is measured in the following process. Glucose as an 
analyte in a blood sample reacts with glucoseoxidase (GOD) to generate a GOD • H2 complex, from which electrons 
are liberated by an electron carrier such as potassium ferricyanide. The free electrons are captured in an electrode to 
calculate a current value, and thus, a glucose concentration is obtained by calculation from a calibration curve, a cali- 
bration formula or the like that is produced previously. 

20 [0004] In an electrochemical measurement, a disposable device called a biosensor is commonly used. A biosensor 
has a reaction layer and an electrode system. 

[0005] A reaction layer comprises a reaction reagent, an enzyme or a matrix to react specifically with an analyte in 
a sample, and provides a reaction site. 

[0006] An electrode system comprises a working electrode and a counter electrode, and by applying voltage to per- 
25 form an oxidation-reduction reaction, the electrode system captures electrons that are generated by the chemical reac- 
tion occurring in the reaction layer, where the electrons are electric signals of a current from an electron carrier to the 
electrode. 

[0007] A biosensor is used in a combination with a measurement apparatus having various functions such as pro- 
viding a certain voltage in a predetermined period, measuring electric signals transferred from the biosensor, and con- 

30 verting the electric signals into a concentration of the analyte. Such a system is called a biosensor system. 

[0008] Among various methods to apply voltage to an electrode system, a method to apply voltage like a rectangu- 
lar wave with respect to time is called a potential step method. In a typical biosensor system, sample feeding to the bio- 
sensor is detected and subsequently, an open circuit or a voltage with substantially no current flow is provided to 
promote a chemical reaction. Alter a predetermined period of time, a fixed voltage is applied to deliver electrons 

35 between the electron carrier and the electrode, i.e., to perform an oxidation-reduction reaction. A state where the bio- 
sensor is provided with a voltage to perform a desired oxidation-reduction reaction is called an excitation state. In gen- 
eral, a current value at an arbitrarily-predetermined point of time during the excitation state is measured, and the current 
value is converted into a concentration of the analyte, using a calibration curve and a calibration formula that are pro- 
duced previously. 

40 [0009] Methods in which plural excitations and plural current measurements are performed are disclosed, for exam- 
ple, in JP-2651278, JP-A-8-304340, and JP-A-10-10130. 

[001 0] JP-265 1 278 provides a method to determine whether a current flowing in a reaction site changes in accord- 
ance with a relationship with a certain Cottrell current 

[001 1 ] JP-A-8-304340 suggests reduction of measurement errors caused by a reduction type intermediate product 
45 that is generated during a storage of the sensor. 

[0012] JP-A-10-10130 provides a method for discriminating a whole blood sample and an aqueous solution as a 
control in a biosensor system. 

[001 3] Measurement results obtainable by using the biosensor can include errors due to various factors. One of the 
factors causing such errors is a sample physical property. 

so [0014] For example, when the sample is whole blood and the analyte is glucose in the blood, the hematocrit value 
(Hct) as a volume ratio of erythrocyte to the whole blood is found to cause errors in the measurement result, and differ- 
ences between individuals are great. The reason appears to be a rising viscosity of the sample. 
[0015] Concentrations of neutral fat and protein in the blood also affect the measurement result. In a measurement 
of an analyte in a blood, sample physical properties, such as blood cells, lipid and protein, will cause measurement 

55 errors. 

[0016] Generally in a conventional method to avoid any influences by such errors, the composition of a chemical 
reaction layer or of an electrode system in a biosensor is improved. For example, JP-A-62-64940 discloses a biosensor 
in which an enzyme to detoxify a measurement-interfering substance is immobilized on an enzyme-immobilizing mem- 
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brane. JP-A-61-3048 discloses a biosensor comprising not only a biocatalytic electrode but an electrode to detect the 
quantity of measurement-interfering substances. JP-A-60-21 1350 discloses a biosensor comprising two electrode sys- 
tems including an electrode system containing an enzyme and a pigment and also an electrode system provided with 
a porous material containing a pigment only However, these biosensors with complicated structures require compli- 
5 cated manufacturing processes, and the manufacturing cost is also raised. 

[001 7] Sensor sensitivity can be another factor causing measurement errors. In many cases, sensor sensitivity var- 
ies from one manufacturing lot of biosensors to another. Manufactures make sensor sensitivity correction tips for the 
respective lots, and ship the tips with their biosensors. 

[001 8] When a lot for a biosensor is changed, a user should correct sensitivity by using a correction tip correspond- 
w ing to the changed lot before he performs an ordinary measurement (JP-A-4-357452). 

[0019] However, such an operation will impose extra work on the user. Moreover, since the range corrected with 
such a correction tip is determined based on the sensitivity during manufacture of the sensor, changes of the sensor 
sensitivity over time, which occur after shipping, will not be corrected. 

[0020] Since an electrochemical measurement includes a chemical reaction, it is affected also by an environmental 
15 temperature and a sample temperature. Namely, an environmental temperature and sample temperature also can be 
factors of measurement errors. JP-2748196 provides a method of correcting nonlinear temperature dependency of a 
chemical sensor. 

[0021] Therefore, the purpose of the present invention is to provide an electrochemical measurement method that 
can decrease work for users and manufacturers without requiring any complicated structure or process of manufacture 
20 ing a biosensor and a measurement apparatus. Such a method can provide highly accurate results by correcting errors 
of concentration of an analyte in a sample. 

DISCLOSURE OF THE INVENTION 

25 [0022] To achieve the purpose, the present invention provides an electrochemical measurement to measure a con- 
centration of an analyte in a sample by using a biosensor having an electrode system and a chemical reaction layer. 
The measurement comprises calculation to obtain as parameters a current value obtained as a result of application of 
a fixed voltage to the sensor after feeding the sample and the ratio of the current value, and calculation of the concen- 
tration of the analyte by using a statistical technique. The parameters are set adequately for each error factor affecting 

30 the measurement results, and error factors can be corrected directly by selecting parameters that will be affected 
greatly by the sample physical properties, or by selecting current values that will be affected greatly by the sensor sen- 
sitivity. 

[0023] For the above-mentioned parameters, for example, the following parameters P1 and P2 are preferably used, 
and the parameters are obtained by applying a predetermined voltage twice to a biosensor after feeding a sample in 
35 order to promote an electrochemical reaction. 

P1 : a ratio (If/lb), where (If) is the value of maximum current or a current occurring after the maximum in a first exci- 
tation, and (lb) is the value of current (lb) read at any point in a second excitation 
P2: a current value (lb) read at any point in the second excitation 
40 With a statistical technique using these two parameters errors caused by sample physical properties and sen- 

sor sensitivity can be corrected with high accuracy. Therefore, the present invention provides measurement with 
high reliability without requiring any complicated structures or methods of manufacturing sensors, or imposing extra 
work on users. 

In addition to the above-mentioned parameters P1 and P2, the following parameters P3 and P4 also can be 
45 used. 

P3: a value (l/AI(y)) obtained by normalizing a derivative value or a differential value of a current at any point in a 
second excitation, with a current value at the same point. 

Here, "normalization" means adaptation of a ratio of either the derivative value or the differential value to a cur- 
rent value at the same point of time in order to make the derivative value or the differential value a parameter not 
so dependent on the concentration of the analyte. When a ratio of the differential value to a current value is adapted, 
the current value can be read at any point between two points having a finite difference. 
P4: a ratio (lb(a)/lb(P)) of an initial current value (lb(a)) to a terminal current value (lb(P)) in a second excitation 

[0024] The inventors consider that these parameters are indices as shown below. However, as it still remains in the 
55 realm of speculation, one parameter will not always correspond to one index. 

P1 : mainly, a current value affected greatry by sample physical properties 

P2: mainly, a current value to indicate a concentration of an analyte in a sample 
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P3: mainly, a current value indicating diffusion of a compound and a mixture of a sample and a reactant inside a 
sensor at a chemical reaction part 

P4: mainly, a current value affected greatly by the electrode sensitivity of a sensor 

5 [0025] In the present invention, the current value lb at any point in the second excitation, which is used for calculat- 
ing the parameters P1 and P2, is preferably the terminal current value (lb(P» in the second excitation. 
[0026J When a value of at least either the parameter P3 or P4 is out of the expectation range, the value is preferably 
substituted by a boundary value of a closest expectation range. This correction will be regarded as an "off-value correc- 
tion". An expectation range in the present invention is a parameter range that is expected from the concentration of an 

10 analyte. 

[0027] It is preferable in the present invention that correction formulas including a plurality of statistical techniques 
corresponding to an environmental temperature or a sample temperature are prepared and an optimum correction for- 
mula is selected therefrom in order to perform an adequate correction corresponding to either the environmental tem- 
perature or the sample temperature. The reason is that the current value obtained by applying voltage to the sensor to 
15 make an excitation state varies considerably by the environmental temperature or the sample temperature, and the 
parameters also vary. 

[0028] Also, when the environmental temperature or the sample temperature is within the boundary region for a 
temperature used in the correction, it is preferable to select both the correction methods adjacent to the boundary in 
order to calculate the respective correction ranges, so that a concentration of the analyte in the sample is calculated 

20 based on a value obtainable by adding either an average value or a weighted average value to a current value to be 
converted into a concentration. When the environmental temperature or the sample temperature corresponds to the 
boundary region, measurement accuracy will be improved further by selecting both the correction methods adjacent to 
the boundary in order to obtain an average value or a weighted average value of the obtained correction value rather 
than selecting one of the methods. 

25 [0029] In the present invention, for example, a discriminant function or a Mahalanobis' distance can be used as a 
statistical technique. 

[0030] The following explanation concerns a case to use a discriminant function as a statistical technique. 
[0031] A discriminant function means a linear expression Z=f(x 1 , x 2 , -..) defined by plural parameters (x 1t x 2 , ...) 
in order to discriminate to which of two previously-produced groups (e.g., G1 and G2) the measurement data belong. 
30 [0032] Z value code is used to discriminate to which of the groups the data belong. For example, a function to dis- 
criminate the level of Hct of a blood sample is defined by using a low Hct group and a high Hct group G 2 as popula- 
tions. 

[0033] It is preferable to obtain a measurement result by preparing plural discriminant functions in accordance with 
a concentration of an analyte in a sample, selecting an adequate discriminant function from the discriminant functions 
35 using lb of the parameter P2 as an index, calculating a discriminant score Z value from the selected discriminant func- 
tion, and by correcting errors with the Z value. 

[0034] When plural discriminant functions are selected by the indices, it is preferable that Z value is calculated from 
each of the discriminant functions, and a correction range corresponding to the Z value is calculated, so that the meas- 
urement result is obtained from either an average value or a weighted average value of the correction range. As men- 
40 tioned above, the measurement accuracy will be improved by using an average value or a weighted average value of 
the correction range rather than using a correction range of any one of the functions. 

[0035] In a correction using the discriminant function, the correction can be performed to determine to which pop- 
ulation the data belong, by using the Z value code. However, it is preferable that the correction range is calculated based 
on the discriminant score Z value. For example, the discriminant score Z value can be divided into a range to be cor- 

45 rected and a range not to be corrected (correction range is 0). 

[0036] The range to be corrected can be divided further into a range to be corrected by a correction range propor- 
tional to the discriminant score Z value, and a range to be corrected by a fixed correction range without reference to the 
discriminant score Z value. This can achieve a correction with higher accuracy when compared to a technique perform- 
ing correction by a grouping depending on the code of the discriminant score. This correction technique will be regarded 

so as "nonlinear correction method" in this description. 

[0037] In the measurement of the present invention, similarly, a measurement with high accuracy can be achieved 
by a selection of Mahalanobis' space and a Mahalanobis* distance D in place of the selection of the above-mentioned 
discriminant function and calculation of the discriminant score Z value as a statistical technique. 
[0038] In the present invention, preferably, errors caused by sensor sensitivity are corrected after errors caused by 

55 sample physical properties are corrected. Since lb as the parameter P2 is considered to reflect the concentration of the 
analyte in the sample, it is preferable that errors caused by the sample physical properties and sensor sensitivity with 
respect to the lb are corrected, and subsequently, the concentration of the analyte in the sample is calculated based on 
a value obtained by multiplying this correction value by a temperature correction coefficient. 
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[0039] Furthermore, a measurement apparatus of the present invention comprises a means to measure either an 
environmental temperature or a sample temperature, a means to detect sample feeding, a means to apply a certain 
voltage at a predetermined point of time, a means to measure a current value generated by an electrochemical reaction, 
and a means to convert the measured value of the current into a concentration of an analyte in a sample, wherein the 
5 temperature measuring means, the means for detecting sample feeding, the applying means, the means for measuring 
current value and the conversion means are controlled to perform the present invention. The control is performed typi- 
cally with a microcomputer having the programmed procedures. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 

[0040] 

F|G 1 illustrates a potential step in an embodiment of the present invention. FIG 2 is a flow chart to illustrate meas- 
urement and arithmetic execution in an embodiment of the present invention. FIG 3 is a detailed flow chart with 
is respect to correction I in an embodiment of the present invention. 

BEST MODE FOR CARRYING OUT THE INVENTION 

[0041] A typical procedure for measurement in the present invention will be explained below with a reference to an 
20 example in which the analyte is glucose in blood, and the explanation is based on FIG 1 illustrating a potential step, FIG 
2 illustrating a measurement-arithmetic execution flow, and FIG 3 illustrating a correction flow. In this example, errors in 
the terminal current value in a second excitation, which are caused by sample physical properties and sensor sensitiv- 
ity, are corrected with a discriminant function, and subsequently, temperature correction is performed and the value is 
converted into a glucose concentration. 
25 [0042] FIG t illustrates a potential step used in the technique of the present invention. A predetermined voltage is 
applied to an electrode system of a sensor after attaching the sensor to a measurement apparatus. Current begins to 
flow when a sample is fed to the sensor. A point of time at which a current of at least a certain level is detected is deter- 
mined as a time of sample feeding, and it is defined as zero second. 

[0043] A fixed voltage is applied after the detection of the sample feeding, and the application is cancelled after t 1 
30 seconds. Therefore, an interval between 0-t 1 seconds is a first excitation. Either an open circuit of at most 3 seconds or 
a process to apply voltage not to cause a substantial oxidation-reduction reaction can be included just after the detec- 
tion of the sample feeding, because in the present invention, what is important is a value of maximum current or a value 
of a current thereafter in the first excitation. Voltage in the first excitation is typically in a range from 300 to 900m V, and 
tj is in a range from 3 to 15 seconds in general. As shown in the figure, the current value in the first excitation initially is 
35 raised and then lowered after having its maximum value. 

[0044] After the first excitation, either an open circuit or a voltage not to cause a substantial oxidation-reduction 
reaction is provided in the interval of t r t 2 seconds. During the interval, glucose, GOD and an electron carrier (e.g., 
potassium ferricyanide) react in a chemical reaction layer of the biosensor. The interval of t-, to t 2 is 3 to 60 seconds in 
general. During t 2 -t 3 seconds, a voltage to cause a desired oxidation-reduction reaction is applied again to the electrode 
40 in order to perform a second excitation. As a result, electrons are delivered between the electron carrier and the elec- 
trode, and a current is observed. A voltage in the second excitation is typically 300 to 900mV, and the time is typically 
2 to 1 0 seconds. The following is a series of the above-mentioned reactions. 

C 6 H 12 O 6 +GOD->C 6 H 10 O 6 +GOD «H 2 

45 

GOD • H 2 +2Fe(CN) 6 3 --*GOD+2Fe(CN) 6 4 >2H + 
FetCNfe^FefCN^+e- 

50 [0045] The above-described four parameters are calculated using a current value measured in the potential step. 
First, a current value (lf(a)) of the maximum or a current thereafter in the first excitation and a current value (lb(P)) in 
the second excitation are measured to calculate P1 (1f(a)/lb(P)) and P2 (lb(P)). The If (a) is obtained typically by meas- 
uring a current value within 15 seconds after the maximum value is observed. The lb(P) preferably is obtained from a 
terminal current value 2-10 seconds after the start of the second excitation (y, and more preferably, it is obtained from 

55 the current value at the end point (ts). 

[0046] P3 (l/AI), which is a value obtained by normalizing either a derivative value or a differential value of a current 
at any point in the second excitation by the above-mentioned current value, is calculated from measured values com- 
prising a current value (lb(y)) at a point (y) in the second excitation and also a current value (lb(y+A)) at a point (y+A) 
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(1) When the Z value is positioned toward the extreme plus side, it is corrected toward the minus side with a fixed 
range. 

(2) When the Z value is positioned toward the plus side, it is linear-corrected toward the minus side. 

(3) When the Z value is positioned in the vicinity of zero, no correction is performed. 

5 (4) When the Z value is positioned toward the minus side, it is linear-corrected toward the plus side. 

(5) When the Z value is positioned toward the extreme minus side, it is corrected toward the plus side with a fixed 
range. 

[0059] When errors caused by sample physical properties is not corrected (correction range A=0), errors caused 
w by sensor sensitivity are corrected. Discriminant functions for correcting sensor sensitivity also are prepared for the 
respective concentration ranges of analyte. The discriminant functions for correcting sensor sensitivity are selected with 
an index of a value obtained by multiplying by lb(P). When a single discriminant function for correcting sensor sensitivity 
is selected, a discriminant score Z value is calculated from the parameter and the discriminant function, and a correc- 
tion range corresponding thereto is calculated. The discriminant function for correcting sensor sensitivity is produced 
15 by using the parameter, as a discriminant function to discriminate populations classified based on differences in the 
electrode sensitivity or the like. When plural discriminant functions for correcting sensor sensitivity are selected, discri- 
minant score Z values and correction ranges corresponding to the Z values are calculated from the respective discrimi- 
. nant functions for correcting sensor sensitivity, and the average value of the correction ranges is determined to be a 
correction range B. 

20 [0060] When the above-mentioned discriminant function (s) for correcting sensor sensitivity is/are used, a correction 
range intermediate value C is determined in the following manner. 

(1) When correction range A*0, 

correction range intermediate value C = correction range A 
25 (2) When correction range A=0, 

correction range intermediate value C = correction range B 

[0061 ] When plural discriminant functions for correcting sample physical properties are selected, a correction range 
intermediate value C is calculated in the above-mentioned manner from the respective discriminant functions for cor- 

30 recting sample physical properties, and the average value or a weighted average value is determined as a final correc- 
tion range D. When a single discriminant function for correcting sample physical properties is selected, the correction 
range intermediate value C is the final correction range D. A value obtained by adding the final correction range D to 
the terminal current value (lb(P)) in the second excitation is the final current value lb(P) tower-middle >in the lower-and- 
middle-temperature side correction. 

35 [0062] In a high-temperature side correction, procedures such as data-taking, production of populations and discri- 
minant functions and calculation of correction range are performed similar to the tower-and-middle-temperature side 
correction, so that a high-temperature side final current value lb(P) highJis calculated. 

[0063] From the temperatures measured with a thermometer, a pre-temperature correction current value lb(p)' is 
obtained. For example, lb(P)'=lb(P) tower-middle > at a temperature from 10 to 30°C, and lb(p)'=lb(P) high ) at a tern- 

40 perature from 35 to 40°C. The temperature range from 30 to 35°C is determined as a boundary region, and either the 
average value or the weight average value of lb(P) tower-middle) and lb(P)high)is determined as lb(p)'. The boundary 
region can be selected experientially by taking the types of the analytes and some other factors into consideration. 
[0064] The lb(P)' is multiplied by a temperature correction coefficient (T1 ) that is calculated previously from an envi- 
ronmental temperature or a sample temperature. The thus obtained value is determined as a post-correction final cur- 

45 rent value Ib(P)" of lb(P). And the lb(P)" is converted into a glucose concentration, using a calibration curve or a 
calibration formula that is produced previously. 

[0065] In the respective sections to select discriminant functions for correction as shown in FIG. 3, only © (©' or 
(3)" for a discriminant function for correcting sensor sensitivity) is selected if a single type should be selected. The cor- 
rection ranges are calculated respectively from ® and © (®' and or <D M and @ M for discriminant functions for 
so correcting sensor sensitivity) if two types are selected, and the average values are determined as the correction range 
intermediate value and the final correction range. 

[0066] The present invention will be explained in detail with respect to specific embodiments. 
[0067] A glucose, sensor system to quantify a glucose concentration in a blood sample will be taken below as an 
example of a biosensor system. 
55 [0068] A glucose sensor used in the glucose sensor system is configured as follows. On an insulating substrate 
comprising polyethylene terephthalate (PET), a carbon electrode system and an electrical insulating layer are formed 
by screen printing, and the carbon electrode system comprises a silver lead, a working electrode and a counter elec- 
trode. The electrical insulating layer has a predetermined area in which the working electrode and the counter electrode 
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are exposed, and partially covers the silver lead. On the thus formed electrode system, a layer of carboxymethylcellu- 
lose (CMC) as a hydrophilic polymer and further, an (enzyme+mediator+CMC) layer are formed. The (enzyme+media- 
tor+CMC) layer comprises glucose oxidase (GOD) as an enzyme, potassium ferricyanide as a mediator (electron 
carrier) and a CMC. The combination of the CMC layer and the (enzyme+mediator+CMC) layer is regarded as a chem- 
ical reaction layer. An insert comprising a cover and a spacer further is formed, so that a fixed amount (about 3uL) of a 
blood sample contacting with the insert is fed by a capillary phenomenon to the chemical reaction layer and to the elec- 
trode system. For the measurement apparatus, a commonly-used potentiostat 100B/W (supplied by BAS) was used. 
While a typical biosensor system uses a combination of a single test cell and an exclusive apparatus, a generalized 
apparatus is used in the present embodiments. However, the embodiments can be performed easily even if a simple 
apparatus of a conventional technique is used. The measurement in the embodiments was performed in an environ- 
mental test lab, but the embodiments can be performed also by using a simple apparatus containing a temperature sen- 
sor. 

[0069] The potential steps were provided in the following manner. 

[0070] A glucose sensor was attached to an apparatus, and a voltage of 500mV was applied to the glucose sensor 
before sample feeding and for 7 seconds after the detection of sample feeding; After an open circuit of 18 seconds, a 
voltage of 500mV was applied again for 5 seconds. The voltage was applied each time for the following purposes. The 
application before the sample feeding is to detect a current generated electrochemically by the sample feeding. Sample 
feeding is detected as a current of at least a certain level, and thus, t=0 will be defined. The 7-second application after 
the sample feeding is for the first excitation, and the main purpose is a detection of difference in the current profile 
caused by the difference in the sample physical properties. The open circuit of 18 seconds is to promote the chemical 
reaction and to store Fe(CN) 6 4 ' as a matrix for an electrochemical reaction. The subsequent 5-second application of 
voltage is for the second excitation, and the main purpose is to convert the Fe(CN)6 4 - to FetCN^electrochemically, 
and to detect a difference in the current profile caused by the difference in the sensor sensitivity and the concentration 
of the analyte's concentration. 

(First Embodiment) 

[0071] First, errors caused by Hct were corrected, since sample physical properties are the most influential factor 
to cause errors in a measurement result. 

[0072] The current was measured at every 0.1 seconds just after the detection of the sample feeding and also after 
the re-application, and current values at the following three points were used for calculation: the peak value l p that was 
measured just after the detection of the sample feeding; 1(0.5) that was measured 0.5 seconds after the re-application 
of current; and 1(5) that was measured 5 seconds after the re-application of current The l p was calculated in the foHow- 
ing manner. Namely, L was l(t) at the moment of l(t)>l(t+0.1) in a current just after the detection of the sample feeding. 
The respective current values were measured for the following purposes. The l p was measured to detect a difference 
in the current profile caused by a difference in the sample physical properties such as a blood viscosity. The 1(0.5) was 
measured to see the ratio to L, so that the glucose concentration factor of l p would be reduced. The 1(5) was measured 
to obtain a current value that depends on the glucose concentration and that can be used for a glucose concentration 
conversion. . 
40 [0073] l(0.5)/l p was obtained for an index to indicate the blood viscosity, and it was determined as a,. In this embod- 
iment, the following experiment was performed indicating erythrocyte amount with hematocrit, since erythrocyte is a 
blood cell in a blood sample and it can be adjusted easily. 

[0074] A relation between a, and 1(5) is shown below, and the relation distinguishes a specimen of 30 /o hematocrit 
and a specimen of 60% hematocrit regardless of the glucose concentration in the blood sample. 

45 

z=-23.1 15xa !+1 .254xl(5)+41 .580 

[0075] The level of hematocrit in the blood sample is reflected by the value of z. When z is positively large, the 
hematocrit in the blood sample is considered to be small. . 
so [0076] A current value l'(5) after the correction of the sample physical properties was calculated by the following 
expression from the relation between and l(5). 

l'(5)=slopex(-z)+l(5) 

55 In this equation, slope=0.0238xl(5)-0.0338. 

[0077] The thus calculated current value f(5) with a corrected blood viscosity will be converted into a glucose con- 
centration value by using a certain calibration curve table. 

[0078] The following data were obtained using a blood sample (glucose viscosity: 536mg/dL, hematocrit: 30%). 
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l p =10.210(nA) 
l(0.5)=20.290(uA) 

5 l(5)=10.110(uA) 

[0079] Regarding the data, a current value P(5) with a corrected blood viscosity can be obtained through the follow- 
ing arithmetic executions. 

10 a l =l(0.5)/l p =1.987 

z=-23. 1 1 5xa 1 +1 .254x|(5)+4 1 .58=8.328 

slope=0.0238xl{5)-0.0338=0.207 

15 

l'(5)=slopex(-z)+l(5)=8.387 

This l'(5) will be converted into a glucose concentration value by using a calibration curve that has been produced. 
[0080] A blood sample used in this embodiment was prepared in the following manner. Blood samples collected 
20 from four persons with the identical blood type were mixed to prepare blood samples having hematocrit of 30%, 45% 
and 60%, respectively. An aqueous solution of glucose was added to these blood samples in order to adjust the glucose 
concentrations to 83mg/dL, 307mg/dL, and 536mg/dL. Measurement was performed using these 9 types of blood sam- 
ples, where the measurement number n was determined to be 1 5. 

[0081] As explained later in detail, errors caused by hematocrit with regard to blood samples were corrected sue- 
25 cessfully in this embodiment. 

[0082] A coefficient of variation (C.V. (%)) is obtained by a calculation based on a group of measurement data for 
blood samples having different hematocrits for the respective glucose concentrations. The following Table 1 shows the 
coefficient of variation to compare the values before and after the correction of the sample physical properties. 

30 

Table 1 . 



35 



Glucose concentration 


Pre-correction 


Post-correction 


83mg/dL 


5.40% 


3.92% 


307mg/dL 


10.69% 


5.67% 


536mg/dL 


16.52% 


8.56% 



[0083] The hematocrit differs depending on individuals. When the standard value is 45% in the present invention, 
40 the divergence (%) of the hematocrit 45% specimens from the measured values for the respective glucose concentra- 
tions are as shown in Tables 2 and 3. Here, divergence of the post-correction data means a divergence from a corrected 
average value of hematocrit 45%. Table 2 shows the divergences (%) from measured average values with respect to 
the cases of hematocrit 30% and of hematocrit 45%. 

45 

Table 2 



50 



Glucose concentration 


Pre-correction 


Post-correction 


83mg/dL 


4.27% 


-0.15% 


307mg/dL 


14.43% 


1.50% 


536mg/dL 


15.69% 


-1.92% 



[0084] It is clear that errors caused by the physical properties of the blood samples were improved as a result of the 
55 correction. 

[0085] Table 3 shows divergence (%) of measured average values for hematocrit 60% from those for hematocrit 
45%. 
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Table 3 
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Glucose concentration 


Pre-correction 


Post-correction 


83mg/dL 


-4.31% 


»1.77% 


307mg/dL 


-10.85% 


-4.10% 


536mg/dL 


-22.21% 


-9.12% 



15 



(Second Embodiment) 

[0087] In this embodiment, a correction of errors caused by an influence of Hrt *q c *m«i a ^ ■ , 

second excitation a ^Tbf^a was Tf™ *? meaSUred " SeC0 " dS *° 

measured 5 seconds after the LI" " "* and 3 «*" lb < 5 > *- «- 

[0089] The following four parameters were used for the operation: 

& P1:lf(7)/lb(5) 
P2:lb(5) 

P3: l/AI(0.5)=lb(0.5)/{lb(0.5)-lb(1 )} 
P4:lb(0.3)/lb(5). 
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II indicates a correction under a high temperature condition. The main reason for classifying the correction method into 
two categories depending on the temperature is that lf(7)/lb(5) is considerably temperature dependent. 
[0093J Four types of discriminant functions for Hct correction with regard to the correction I were produced: func- 
tions for 30-75mg/dL, 75-1 50mg/dL, 1 50-300mg/dL, and 300-500mg/dL In a case of 75-1 50mg/dL, a discriminant func- 

5 tion to discriminate two populations of H t and Lj for 75mg/dL and 1 50mg/dL was obtained. In this case, the discriminant 
function was produced so that the discriminant score Z value of a data having high values due to the influence of Hct 
(i.e., a data of high Hct specimens) would be plus. In producing the populations, data including parameters regarded as 
abnormal values as a result of a null hypothesis test of (Grubs-Smirnor were eliminated, since the populations were 
preferably as uniform as possible. 

w [0094J The discriminant function for Hct correction was produced by selecting a value as an index obtained by mul- 
tiplying a temperature correction coefficient (T1) by lb(5), where the temperature correction coefficient was calculated 
previously from the environmental temperature, and the selection numbers were 1 -2. f 1 was determined as 1 at 25°C, 
a temperature that is considered as a standard. A case where plural discriminant functions were selected will be 
described later. When a single discriminant function for Hct correction was selected, the parameter was substituted for 

75 the discriminant function, and a discriminant score Z value and a correction range A corresponding thereto were calcu- 
lated. The correction range A was determined with a nonlinear correction method. Namely, when 
z a <z b <z c <z d(Za.<Zb.<Zc,<Z d are constants experientiatly determined for every discriminant function) and when a 
current value lb(5) and a discriminant score Z value were obtained, corrections were performed as follows: when Z<Z A , 
correction was performed with a fixed range toward the plus side regardless of Z value; when Z A <Z<Z B , a linear correc- 

20 tion was performed toward the plus side corresponding to the Z value; when Z B <Z<Z C , no correction was performed; 
when Z C <Z<Z D , a linear correction was performed toward the minus side corresponding to the Z value; and when Z D <Z, 
a correction with a fixed range was performed toward the minus side regardless of the Z value. The correction with a 
fixed range means correcting with a fixed range determined by a measured lb(5) for every discriminant function. A linear 
correction is determining a correction range from a proportional relation with the Z value with an upper limit of the 

25 above-mentioned fixed range. For example, if a linear correction is performed toward the minus side when 5<Z<1 0, the 
correction range is -0.10 for Z>10, a correction range for Z=8 is calculated according to the following equation (1). 

Correction range A=-0.10x(8-5)/(10-5)=-0.06 (1) 

30 [0095] The term 'correction of a fixed range' means that a correction range will have a fixed value once lb(5) is 
determined. The correction range itself varies in accordance with lb(5). When no Hct correction was performed (correc- 
tion range A=0), the sensor sensitivity as mentioned above was corrected as described below in order to obtain a cor- 
rection range B. When Hct correction was performed (correction range A*0), correction of sensor sensitjvity will not be 
performed. 

35 [0096] Seven types of discriminant functions for correcting sensor sensitivity with regard to correction I were pre- 
pared: functions for 75-1 OOmg/dL, 100-150mg/dL, 150-200mg/dL, 200-250mg/dL, 250-300mg/dL, 300-400mg/dL and 
400-500mg/dL. In a case of 75-1 OOmg/dL, a discriminant function to discriminate two populations of H 2 and L 2 for 
75mg/dL and 1 0Omg/dL was obtained and the function was determined as a discriminant function for correcting sensor 
sensitivity for the above-identified concentration ranges. In this case, the discriminant function was produced, so that 

40 the discriminant score Z value measured with a highly sensitive sensor would be plus. In producing the populations, 
data including parameters regarded as extremely high/low as a result of a null hypothesis test of (Grubbs-Smirnov were 
eliminated, since the populations were preferably as uniform as possible similar to the above-mentioned case of pro- 
ducing a discriminant function for Hct correction. 

[0097] The discriminant function for correcting sensor sensitivity was selected with lb(5) as an index, and the selec- 
45 tion numbers were 1-2. When a single discriminant function was selected, the parameter value was substituted for the 
discriminant function, and a discriminant score Z value and a correction range corresponding thereto were calculated. 
The correction range was regarded as a correction range B. The correction range was determined with a nonlinear cor- 
rection method similar to the case of Hct correction. When plural discriminant functions were selected, correction 
ranges were calculated from the respective discriminant functions, and the average value was regarded as the correc- 
so tion range B. 

[0098] A correction range intermediate value C was determined as follows when a single discriminant function for 
Hct correction was selected. 

(1) When correction range A*0, correction range intermediate value C=correction range A 
55 (2) When correction range A=0, correction range intermediate value C=correction range B 

[0099] When two types of discriminant functions for Hct correction were selected, a correction I final correction 
range D was obtained in a similar way by obtaining an average value or a weighted average value of the correction 
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range intermediate values C obtainable from the respective discriminant functions for Hct correctioru When only one 
type of discriminant function for Hct correction was selected, a correction range intermediate value C thereof was deter- 
mined as the correction I final correction range D. The correction I final current value lb(5), was a value obta.ned by add- 
ing the final correction range to lb(5). 

[01001 Regarding the correction II, the technique of the correction I was adapted to take data, produce popu ations 
and discriminant functions, and to calculate correction ranges, so that a Correction II final current value lb(5)„ was 

[OW] 6 A pre-temperature correction current value lb(5)' was determined as follows from a temperature T (*C) 
detected with a temperature sensor. 

WhenT<30, lb(5)'=lb(5) , . 
When 30<r<35, lb(5)'=(35-T)/(35-30)xlb(5) ,+CT-30)/(35-50)xlb(5) „ . 
WhenT>35, lb(5)*=lb(5) „ . 

15 [01 02] For example, lb(5)' can be obtained from the following equation when T=32°C: 

lb(5) 4 =3/5xlb(5) ,+2/5xlb(5) , t 

r0103] lb(5) M was obtained by multiplying thus obtained lb(5)' by a temperature correction coefficient (T2)^ This 
20 lb(5)" was converted into a glucose concentration value according to a calibration curve table that was produced sep- 

P104]* For example, the following data were obtained by measuring blood specimens of a glucose concentration of 
1 54mg/dL and Hct 65% under an environmental temperature of 25°C. 

25 lf(7)=2.142(nA) 

!b(0.3)=7.578(iiA) 

lb(0.5)=6.312(nA) 

30 . v 

lb(1)=4.932(nA) 

ib(5)=2.142(uA) 

35 [0105] The following parameters are obtained from the above-identified data. 

If(7)/lb(5)=0.8069 
l/AI(0.5)=lb(0.5)+{lb(0.5)-lb(1 )}=4.575 

40 

lb(0.3)/lb(5)=2.855 
lb(5)=2.654 

45 [0106] At this time, expectation ranges of l/AI<0.5) and lb(0.3)/lb(5) are calculated from the value of lb(5) ^respec- 
tively to be 3.050-4.496 and 2.817-4.361. Since l/Al(0.5) exceeds the upper limrt of the expectation range and thus is 
considered as an off-value, l/AI(0.5)=4.496 in the subsequent calculations. 

[0107] Correction I is performed. A discriminant function for Hct correction is selected with an index of lb(5)-2 .65* 
When lb(5) ranges from 2.407 to 3.400, a discriminant function for 75-150mg/dL and a * scnm ^ 
50 300mg/dL are used. The discriminant functions are defined respectively as follows. A discriminant function for 75- 
1 50mg/dL • Hct correction: 

Z(L)=34.525xP1-27.333xP2-21.717xP3+3.1870xP4+141.90. 

55 A discriminant function for 1 50-30mg/dL • Hct correction: 

Z(N)=331.25xP1-16.739xP2-22.331xP3+6.7191xP4-158.00. 
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In these equations, P1=lf(7)/Ib(5) , P2=l/A 1(0.5) , P3=lb(0.3)/lb(5) , and P4=lb(5) . 

[0108] In this embodiment, Z(L)=-5.753. Since a linear correction is performed toward the plus side when - 
8<Z(L)<-4, a correction range A L is calculated from the following equation: 

5 A L ={0. 1 205xib{5)-0. 1 049}x{-4-Z}/{-4-(-8)}=0.1 024. 

[0109] The sensor sensitivity will not be corrected since A^O. The value of 0.1 024 will be a correction range inter- 
mediate value C L . 

[01 10] In this embodiment, Z(N)=-1 1 .887 . Since no correction is performed when -1 5<Z(N)<5, the sensor sensitiv- 
w ity will be corrected. Since lb(5)xtemperature correction coefficient (T1 )=2.654 and this figure is within a range of 2.582- 
3.316, a discriminant function used for correcting sensor sensitivity is of 1 50-200mg/dL. The discriminant function of 
1 50-200mg/dL is defined as follows. A discriminant function for correcting sensor sensitivity of 1 50-200mg/dL: 

Z(S 3 )=43.581 xPI-25.041 xP2-25.472xP3+2.0536xP4+141 .86. 

IS 

In this equation, P1=lf(7)/lb(5) , P2=l/AI(0.5) , P3=lb(0.3)/lb(5) , and P4=!b(5) . 
[01 1 1 ] In this example, Z(S 3 )=-2.827 . 

[0112] Since a linear correction is performed toward the plus side when -6<Z(S 3 )<-2, a correction range B N is cal- 
culated from the following equation. This B N will be a correction range intermediate value C N . 

20 

B N =C N ={0.03 1 2xlb(5)-0.01 21 }x{-2-Z(S3)}/{-2-(-6)}=0.01 46 

[0113] Since two discriminant functions for Hct correction are selected in this example, a final correction range D, 
of the correction I is an average value obtained from the two correction range intermediate values C L and C N : 

25 

D ,=(0.1024+0.0146)^2=0.0585. 

Therefore, a post-correction current value lb(5)j obtainable by the correction I can be obtained in the following manner. 

30 lb(5) ,=2.654+0.0585-2.713 

[0114] Since the temperature is 25°C and the pie-temperature correction current value lb(5)'=lb(5) , , calculation 
procedures for obtaining lb(5)u will be omitted from the following description. Since the temperature is 25°C, the tem- 
perature correction coefficient is 1 and the post-temperature correction current value lb(5) ,, =lb(5)'=2.713 . This lb(5)" 
35 value is substituted for a calibration curve produced separately in order to be converted into a glucose concentration. 
An average value of lb(5) of a data of a glucose concentration of 154mg/dL and Hct 45%, i.e., an expectation value of 
lb(5) for a glucose of 154mg/dl was 2.790. A correction according to the present invention decreased divergence from 
the expectation range (pre-correction value=2.1 42; post-correction value=2.713). 

[0115] In this embodiment, errors caused by sample physical properties and sensor sensitivity were decreased. 
40 This will be explained below referring to an result of a case where the glucose concentration was 1 50mg/dL 

[0116] Ten lotsxnumber N (10) of sensors were measured for every Hct specimen. The coefficient of variation 
(CV)(%) for the entire measurement data varied before and after the correction as described below. The table below 
shows that errors caused by sensor sensitivity is decreased. 

45 

Table 7 



Specimen Hct 


Pre-correction 


Post-correction 


Correction effect (Pre- 
correction-Post-correc- 
tion) 


25% 


4.7% 


4.3% 


A0.4% 


45% 


4.5% 


4.2% 


AO.3% 


65% 


5.1% 


4.7% 


AO.4% 



55 

[0117] Since the hematocrit differs depending on individuals, the standard value is 45% in the present invention. 
Divergences (%) of the respective Hct specimen current average values from a Hct 45% specimen current average 
value are shown in Table 2. It is clear from this table that errors caused by Hct were corrected further. 
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Table 8 
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Specimen Hct 


Pre-correction 


Post-correction 


Correction effect (Pre- 
correction-Post-correc- 
tion) 


25% 


15.6% 


8.2% 


A7.4% 


65% 


-7.2% 


-3.8% 


A3.4% 



[01 18] As an index to indicate a maximum divergence of an obtained measurement value from a true value or from 
an expectation value, a total error (TE) is considered. TE is defined as follows. 

15 TE(%)=DEV(%)+1.65xCV(%) 

[01 1 9] A coefficient 1 .65 of a CV term is a value of Z value (a value not as a discriminant score but as a value in a 
standard normal distribution) composing one-sided probability of 5% of a standard normal distribution, and TE defined 
in this term is a divergence range from a true value to which 95% of the entire data is considered to be included 
20 [01 20J TE in this example is as shown in Table 3. 



Table 9 



25 



30 



Specimen Hct 


Pre-correction 


Post-correction 


Correction effect (Pre- 
correction-Post-correc- 
tion) 


25% 


23.3% 


15.3% 


A8.0% 


45% 


7.4% 


7.0% 


AO.3% 


65% 


15.6% 


11.5% 


A4.4% 



INDUSTRIAL APPLICABILITY 

[01 21 ] As mentioned above, according to the measurement method of the present invention, errors with respect to 
analytes caused by factors such as sample physical properties and sensor sensitivity can be corrected efficiently with- 
out modifying compositions of a chemical reaction layer or an electrode system of a sensor, and thus, a reliable meas- 
urement value can be obtained. 
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Claims 
1 



An electrochemical measurement method to measure a concentration of an analyte in a sample by using a biosen- 
sor, comprising calculation of parameters from a current value that is obtained as a result of application of a fixed 
voltage to the biosensor after sample feeding, and correction of errors of a measured value using these parameters 
by a statistical technique. 

2. The measurement method according to claim 1 , wherein parameters that are affected by errors caused by sample 
physical properties are selected. 

so 3. The measurement method according to claim 1 , wherein parameters that are affected by errors caused by sensor 
sensitivity are selected. 



4. The measurement method according to claim 1 , wherein the concentration of the analyte is obtained by applying a 
predetermined voltage twice to the biosensor after sample feeding, calculating parameters Pi and P2 from the cur- 
rent values provided by the application, and by correcting the errors from these parameters by a statistical tech- 
nique, where 

Pi: a ratio (If/lb), where (If) is the value of maximum current or a current occurring after the maximum in a first 
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excitation, and (lb) is the value of current read at any point in a second excitation; 
P2: a current (lb) read at any point in the second excitation. 

5. . The measurement method according to claim 4, wherein the current value lb at any point in the second excitation 
5 used for the parameters P1 and P2 is a terminal current value (lb(P)) in the second excitation current. 

6. The measurement method according to claim 4, wherein at least one of parameters P3 and P4 is used together 
with the parameters P1 and P2, where 

w P3: a value (l/AI(5)) obtained by normalizing a derivative value or a differential value at any point in the second 

excitation with a current value at the same point; and 

P4: a ratio (lb(a)/lb(p)) of an initial current value (lb(a)) to a terminal current value (lb(P)) in a second excitation. 

7. The measurement method according to daim 6, wherein at least one of the parameters P3 and P4 is substituted 
15 with a closest expectation value when the parameter is out of an expectation range. 

8. The measurement method according to claim 1 , wherein an adequate correction is performed in accordance with 
an environmental temperature or a sample temperature by preparing a plurality of correction formulas using a sta- 
tistical technique corresponding to the environmental temperature or the sample temperature and by selecting an 

20 optimum correction formula from the plural correction formulas. 

9. The measurement method according to claim 8, wherein the concentration of the analyte in the sample is obtained 
when the environmental temperature or the sample temperature corresponds to the boundary region of the tem- 
perature range used for the correction formula, by selecting both the correction formulas adjacent to the boundary, 

25 calculating correction ranges from the respective formulas, adding the average values or weighted average values 
to the current value for concentration conversion, and obtaining the concentration depending on the additional 
value. 

10. The measurement method according to claim 1, wherein a discriminant function is selected from the parameters. 

30 

11. The measurement method according to claim 10, wherein a discriminant score is obtained from the discriminant 
function and the concentration of the analyte is calculated from a correction range proportional to the discriminant 
score. 

35 12. The measurement method according to claim 10, wherein a range to calculate the concentration of the analyte 
depending on a fixed correction range regardless of the discriminant score, and a range to calculate the concen- 
tration of the analyte without performing a correction are set. 

13. A measurement apparatus comprising a meter for measuring an environmental temperature or a sample tempera- 
40 ture, a detector for detecting sample feeding, a device for applying a certain voltage to a biosensor at a predeter- 
mined point of time, a meter for measuring a value of current generated by an electrochemical reaction, and a 
converter for converting the measured value of the current into a concentration of the analyte in the sample; 
wherein the meter for measuring temperature, the detector for detecting sample feeding, the device for application, 
the meter for measuring the current value and the converter are controlled to perform the measurement according 
45 to any of claims 1 to 1 2. 
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